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o\ : 

t^)' ABSTRACT 

^3 : 

<: 

^ . The observed braking index riohs which had been determined for a few young 

pulsars, had been found to differ from the expected value for a rotating magnetic 

^ . dipole model. In addition, the observational jerk parameter, determined for two 

^ . of these pulsars, disagrees with the theoretical prediction niobs = 15 in both 

O . cases. We propose a simple model able to account for these differences, based 

^ . on a growth of the "torque function" K = —17/^2", under the constraint that 

. Uobs is a constant. We show that there is observational evidence supporting the 

^ , latter hypotesis, and derive initial values for several physical quantities for the 

P_|. four pulsars whose Uobs have been measured. 

o ■ 

^ I Subject headings: Pulsars : general — Pulsars : individual : (PSR B0531-I-21, 



X 



PSR B0540-69, PSR B0833-45, PSR B1509-58) 

1. Introduction 

The well-known vacuum dipole model expresses the external torque acting on pulsars 



as 



in = — '^B^R^sin^an'' (1) 
3 

where a is the angle between B and Q, B and R are the magnetic field and the radius of 
the star respectively, I is the star's moment of inertia and c is the velocity of light. The 
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braking index n is usually assumed to be 3 as predicted by the magnetic dipole model 

(see e.g. Manchester & Taylor 1977), but actually the four (young) pulsars with the best 
determinations of braking index show values < 3 (see Blandford & Romani 1988 and 
references therein). In the general literature all these quantities are taken to be constants 
(except of course for VL and its derivatives), which allows one to rewrite eq.(l) as 

h ^ -K^r (2) 

where is a function that absorbs the structural factors, and we will refer to throughout 
this work as the "torque function". 

There is evidence that in 1975 and 1989 the spin rate Q of the Crab pulsar suddenly 
increased by amounts AQ/fi ~ 10~^ and after that continued to spin-down at a faster rate 
(that is, the pulsar continued to spin slower than before) of /STL/VL ~ 10~^ (GuUahorn 
et al. 1977, Lohsen 1981, Lyne, Smith & Pritchard 1992), characterizing a permanent 
deficit in Ct. The same feature is also present in the 1969, 1981 and 1986 glitches (Lyne & 
Pritchard 1987, Lyne, Pritchard & Smith 1993). From cq.(l) it is clear that the peculiar 
events of 1975 and 1989 require either a reduction of / (Alpar & Pines 1993), an increase of 
the magnitude of the magnetic field B (Blandford, Applegate & Hernquist 1983, Muslimov 
& Page 1996, Camilo 1997) or an increase of the angle a (Macy 1974). Other dynamical 
models in which the angle between the magnetic dipole and the rotation axis of the pulsar is 
allowed to vary hace been proposed by Link & Epstein (1997) and Allen & Horvath (1997). 
In our previous work, we assumed several laws for the growth of the angle a and verified 
which ones provided consistent solutions. It was shown that a exponential angular growth 
with a e-folding time of ~ lO'^yr could fit the small braking indexes and jerk parameters 
of the Crab, Vela and B0540-69 pulsars, and although PSR B1509-58 was best-fitted by a 
logarithmic law, an exponential one is not ruled out because of the large uncertainties (25 
%) on its jerk value. 

The key feature here is the variation of the torque function K, which cannot be 
held constant, but instead evolves, growing on timescales shorter than the pulsar's life 
span. This work presents a simple general model for torque function growth, based on the 
observational feature riobs — constant (see also Ruderman 1993 and references therein). 
In the next Sections we will show that this simple hypothesis is enough to determine the 
general behaviour of the pulsar dynamics, along with its evolution on the P x P diagram, 
yet taking into account the low braking indexes displayed by young pulsars. 



2. Consequences of nobs — constant 
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2.1. General model 

To relate the model to the observations, we need to define some physical quantities. 
The first one is the variation of the torque function AK due to persistent shifts /SCt/Ct and 
AVt/Vt, which is easily found from eq.(2) to be 

Dividing eq.(3) by At, here defined as a typical timescale between glitches, we obtain 
a mean increase rate which we shall denote as (■^)- Note that even the exact definition 
of the mean is not completely satisfactory given the different coverages and observational 
biases for each pulsar (we will return to this point below). Using the data from Lyne, 
Pritchard & Smith (1993) (hereafter LPS), it can be found that the mean torque function 
variation caused by glitches {discrete variation) for the Crab is 



with At — 4.6 yr. 

It can be shown that even if the power of ft in the torque expression is exactly 3, the 
observed braking index defined as Uobs — ^^/Cl^ is 

no,. = 3 + 2|, (5) 
and the observed jerk parameter defined as rriobs — Clfl'^/Cl^ is 



rriobs = ^{^riobs - 4:) + {uobs - ^) (^) 

Eq.(5) can be inverted to find the (continuous) torque function variation needed to 
account for the observed braking index of Crab 

§ = (nobs - 3) ^ ~ 1.9 X W-^yr-\ (7) 

So, the discrete contribution is only ~ 15% of the continuous torque function variation. 
That means that even pulsars that have not shown any glitches could have their torque 
functions varying in a continuous way, and also satisfy Uobs < 3. In fact, PSR B1509-58 
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and PSR B0540-69 have never displayed glitches, yet their braking indexes are 2.8 and 2.0 

respectively, although the hypothesis that they are actually glitching at a rate comparable 
to the Crab is not ruled out. Application of eq.(7) to Vela, PSR B1509-58 and PSR 
B0540-69 yields 6.98 x 10-^ 5.25 x 10"^ and 28.8 x 10'^ yr-\ respectively. 

With the expression for Ugbs we derive a form for the jump of the latter in each glitch 
event, namely 



/An Ad AQ\ 
An.,. = n.,,(^— -2_ + _j. (8) 

The division of eq.(8) by At yields a mean variation rate of riobs, which can be 
calculated from LPS for the Crab to be 



(^) ^ -l-5xl0-^r-^ (9) 

assuming that Cl did not vary in the events, as assumed in all data analysis performed 
until now. According to LPS the data show that Uobs appears to be constant within 0.5% 
(~ 0.01) over 20 years, implying \hobs\ < 5 x 10^'^ yr^^, which is consistent with eq.(9). 
Thus we are led to explore a scenario in which pulsars evolve along a constant braking 
index value ^ 3. 

A direct measurement of fiobs in term of observables can be obtained as 



nobs — ^obs ^ + 7T — 2 — = : . (10 

Vfi n nJ n/n ^ ' 

Setting fiobs — 0, we can express rUobs as function of Uobsi ^ and Q, and replace it in 
eq.(6), yielding 



K K Uobs ~ 1 

— ■■ — = 1 — = constant (11) 

which can be now integrated from an arbitrary point to the present values (indicated by 
the subscript p), resulting in 



^-u^y—- (12) 



Eq.(5) can be also integrated in the same way, and we obtain 



Again, integration is performed and with the help of eq.(7) the time dependence of K 
is found to be (for Uobs 7^ 1) 



K 



I + {nobs - l){tp - t)^ 



(14) 



where tp is the true age of the pulsar; the time dependence of Q is trivially recovered 
substituting eq.(14) into eq.(13) 



{riobs - l)(tp - t)^ 



(15) 



It is worth remarking that eq.(15) has the same form of the standard calculation made 
by assuming that the torque function does not change, with Uobs replacing n — 3 (see 
Manchester & Taylor 1977). In other words, it appears as though the actual external torque 
acting on pulsars has the form Cl oc $7"°*". Therefore, several approximations customarily 
made in the literature are exact in the limit Uobs — constant. 



Finally, from eq.(lO) we write 



mobs = nobsC^nobs - !)• (16) 

Again, we note that eq.(16) has the same form as the conventional relation, if we 
replace n = 3 by nobs- 

It is important to note that the lower nobs is, the closer the initial period is to the 
present one. Therefore, the use of the characteristic age r = b^~i) |^ (obtained from 
eq.(15)), though strictly valid when nobs — constant, introduces a non- negligible error 
because Q ~ Qo- 

The important fact is that postulating fiobs = (as suggested by observations) 
completely determines the dynamical evolution of the pulsar, whichever model of variation 
one chooses between magnetic field, moment of inertia or angle a. In Table 1 we show the 
correspondence between K and its derivatives for all specific models. 

In Fig.l we depict the evolutionary tracks of the Crab and B 1509-58 pulsars predicted 
by our model, compared to the standard model tracks, from birth to 5 x 10^ years. Fig.2 
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shows the same for Vela and B0540-69. Remarkably, P for Vela is actually increasing with 
time, although this does not mean a rising torque, because of the compensating effect of 
the power of Q. 

Now, it is easy to calculate the initial period for any pulsar whose Uobs has already 
been determined. In Table 2 we show observed and calculated quantities for the above 

pulsars, where L = IQCI is the energy-loss rate of the rotating magnetic dipole. The data 
used in calculations came from LPS (Crab), Kaspi ct al.(1994) (PSR B1509-58), Taylor et 
al. (1995, unpublished, sec Taylor, Manchester & Lync 1993) (PSR B0540-69), Lync ct 
al. (1996) and Taylor, Manchester & Lync (1993) (Vela). Except in the case of Crab, the 
conventional characteristic age Q/{2Q) was arbitrarily used as the true age of the pulsar, 
lacking a more reliable determination. 

2.2. Crab Pulsar (PSR B0531-|-21) 

The Crab pulsar is known to be among the most active glitching pulsars. If the glitch 
activity has been aproximately the same since its birth, nobs would have decreased by an 
amount {Anobs/ ^t) x ~ 0.15, or 6% of the present value. In the 25 years that have 
passed since the first Crab measurements, ghtches could have provoked a reduction of only 
0.004 (0.2%) in Uobs, far below the observational hmit (see LPS). About 20 years' more data 
will be needed to unambigously detect Ugbs variation. 

The values obtained assuming Ugbs — constant (Table 2) are very close to those 
observed, which is as expected, since Ct was in fact estimated by LPS using eq.(16). The 
good fit to the data reinforces that riobs is zero or very small. 

2.3. Vela Pulsar (PSR B0833-45) 

Some authors (Aschcnbach, Egger & Triimper 1995, Lyne et al. 1996) have discussed 
the possibility of Vela being up to 3 times older than its conventional characteristic age. As 
a test of the sensitivity of this model to the true age adopted, an age doubhng alters the 
values in Table 2, as follows 

Po ~ 26 ms 



Ko ~ 1.56 X 10"'^ s 



(17) 
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Po ~ 60.4 X 10~ 




It is worth to noting that Vela is also one of the most glitch-active pulsars known, 
and, like the Crab, there could be a discrete torque function growth, as discussed in the 
introduction, which would make riohs 7^ 0. However, timing analyses made to date (Lyne 
et al. 1996 being the most recent) have failed to unambigously detect permanent shifts in 
because the typical (non-permanent) shift Afi/Q is 10~^ and the mean time between 
glitches is ~ 2 yr, leading to the conclusion that if such a permanent component does exist, 
it must be negligible. 

3. Conclusions 

We have shown that low braking indices, observed in young pulsars, can be attributed 
to a variation of the torque function K. Besides, there are evidence that Uohs — constant 
in some of these pulsars, and this hypothesis leads to a complete determination of 
their dynamics, which we have calculated. Our expressions are almost identical to the 
conventional ones, when n substitutes nobs, justifying the customarily estimates. 

Prom Table 2, we can see that K must have increased by 25% to 100% since the pulsar 
birth. From the equivalence in Table 1, we verify that this means a reduction of 20-50% in 
the moment of inertia, which can be accomodated neither in the current vortex creep model 
nor as loss of oblateness. So a decrease in the moment of inertia cannot be invoked as the 
main physical agent behind torque evolution, although it probably is in the case of glitches. 

An inspection of Table 2, Fig.l and Fig. 2 reveals sistematically larger initial periods 
and smaller period derivatives than those obtained in the standard model. It becomes clear 
that, within nobs = constant models, the characteristic age is poorly related to the true age 
of pulsars below 10^ yr. This may affect statistical studies of pulsar populations. 

We think this torque evolution, especially if driven by angular growth, probably 
stops at ages of the order of 10^ years, when a = 90" (Allen & Horvath 1997). Further 
improvement will require observations of the braking indices of middle-aged (~ 10^ years) 
pulsars. 
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Table 1. Equivalence between expressions of general and specific models 



General Angle Magnetic Field Moment of Inertia 



K sin^o; B'^ I'^ 

k/K 2a/ tan a 2B / B -I /I 

k/k 2cos(2a)d + (a/d) {B/B) + {B/B) (-2///) - (J//) 

kK/k^ tana[cos(2a) + (ci/2Q;2)] ^{1 + {BB/B^)) 2 + {il/P) 
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Table 2. Comparison between observations and calculated values for young pulsars 





Crab 


PSR B1509-58 


PSR B0540-69 


Vela 




915 


1550 


1660 


11000 


'^obs 


2.5179 ±0.0001 


2.837 ±0.001 


2.04 ±0.02 


1.4 ±0.2 


fiobs [10~^|/r-^] 


-2.7 ± 1.3 


-4± 11 






mobs 


10.23 ±0.03 


14.5 ±3.6 






TTlobs if nobs = 


10.160 ±0.001 


13.26 ±0.01 


6.28 ±0.16 


2.8 ±0.85 


Po [ms] 


19 


39 


25 


52 


Po[lO-''] 


559 


4760 


492 


91 


Lo/Lp{I — const) 


6.7 


180 


8.4 


3.7 


KollQ-^Ks] 


10.8 


186 


12.3 


4.73 


Ko/K, 


0.770 


0.802 


0.510 


0.424 


# {nobs = 0) 


-2.15 


-10.27 


-0.081 


0.737 
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Fig. 1. — In this diagram, the Crab pulsar track is marked with x, while PSR B1509- 
58 is marked with +. Each step is lO^t/r, from birth to 5 xlO^. Squares indicate the 
present situation. Both tracks are intercepted by lines representing the same tracks after the 
conventional model. 
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Fig. 2. — In this diagram, the Vela pulsar track is marked with x, while PSR B0540- 
69 is marked with +. Each step is lO^t/r, from birth to 5 xlO^. Squares indicate the 
present situation. Both tracks are intercepted by lines representing the same tracks after the 
conventional model. 



